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We calculate the entropy produced in the decoherence of a classical field configuration and compare 
it with the entropy of a fully thermalized state with the same energy. We find that decoherence alone 
accounts for a large fraction of the equilibrium entropy when each field mode is only moderately 
occupied. We apply this to theories of relativistic heavy ion collisions, which describe the initial 
state as a collection of coherent color fields. Our results suggest that decoherence may partly explain 
the rapid formation of a high entropy state in these collisions. 



Recent measurements of anisotropic flow in Au+Au 
collisions at the Relativistic Heavy Ion Collider (RHIC) 
point to a very rapid formation of nearly equilibrated 
hot matter in these processes 0, @ . Hydrodynamic cal- 
culations of the collective flow pattern yield estimates 
for the equilibration time r oq < 1 fm/c. This is difficult 
to understand in microscopic theories of thermalization 
0, 0] , which suggest that the thermalization of the dense 
matter proceeds through radiative processes and takes 
several fm/c. 

The initial state of a relativistic heavy ion collision is 
characterized by a highly coherent configuration of quark 
and gluon fields. For the processes contributing to the 
formation of matter near the center of momentum of the 
two nuclei, coherent gluon fields at small values of the 
Bjorken variable x are most important. These fields are 
generated by the quasistatic color charges of the valence 
quarks of the nuclei and can be approximated as ran- 
domly oriented, quasiclassical color fields, often called 
a color glass condensate During the nuclear collision 
gluons are scattered out of this coherent field with a prob- 
ability that is predicted to be close to unity [y, [j| ■ After 
their liberation, the gluons scatter off each other and ra- 
diate additional gluons until they reach an equilibrium 
|. 

Parton cascade models |8| describe the liberation and 
rescattering of gluons and other partons in a probabilistic 
framework based on the relativistic Boltzmann equation, 
starting from an incoherent ensemble of partons. On the 
other hand, models based on an initial state of the color 
glass condensate type treat the equilibration process as 
nonlinear classical evolution of the initial random, but 
coherent color fields 0, @, E3 ■ While entropy production 
in statistical transport theories, such as the Boltzmann 
equation, is a well studied and understood phenomenon, 
the production of entropy by the decoherence of classical 
fields is less well understood. One mechanism for the 
production of entropy is the pair creation of particles in 
strong (chromo-)electric fields |llj . Another mechanism 
is the dynamical chaos generated by the nonlinear field 
equations |l2j . where the Kolmogorov-Sinai (KS) entropy 



describes the apparent rate of entropy production. 

Here we are not concerned with the microscopic de- 
scription of the production of entropy; instead, we ad- 
dress the question of the relative contribution to entropy 
production by (a) the decoherence of an initially coher- 
ent field configuration and (b) the rescattering among in- 
coherent particle-like field excitations, which ultimately 
leads to equilibrium. Since neither the initial color field 
configuration in a fast moving nucleus nor its dynami- 
cal evolution after the collision of two nuclei is very well 
known, we start with a simple case, for which the rel- 
evant calculations can be performed exactly, but which 
is sufficiently general to permit conclusions that can be 
applied to heavy ion reactions. 

The idea that decoherence may play a major part in 
entropy production in heavy ion collisions is not new. 
The notion of a large nucleus acting as a "phase filter" , 
decomposing the quark-gluon wavefunction of a hadronic 
projectile into its incoherent components, was suggested 
over a decade ago 0, Q. The principles of decoher- 
ence of hadronic wavefunctions were investigated exten- 
sively by Elze ^ij in the mid-1990s. The formulation of 
models for the small Bjorken- a; components of hadronic 
wavefunctions as superpositions of classical color fields, e. 
g. the color glass condensate (CGC) model @, now pro- 
vides the theoretical basis for a more concrete treatment. 
In this framework we need to explore the effects of the 
decoherence of quasiclassical, coherent color fields, which 
are present in the nuclei before the onset of a collision. 

The quantum mechanical analogue of a classicical field 
is a coherent state (lif : 



\^[J\) = n eXP ( iakAa kA - *«kA a kA)|0), 



(1) 



where the amplitude akA is determined by the classical 
current J creating the field: 



QfkA = (fiwkV) 1/2 £kA • J(k,07 k ) 



(2) 



Let us begin by considering a single mode kA. The co- 
herent state can be written as a superposition of particle 



2 



number eigenstates: 



|a) = e-HV^ 

n=Q 



(3) 



Being a pure quantum state, \a) is described by a density 
matrix 



(m|a)(Q;|n), 



(4) 



which satisfies the relation p 2 = p and has no entropy: 
S = — Tr/5 hip = 0. 

Complete decoherence of this quantum state corre- 
sponds to the total decay of all off-diagonal matrix ele- 
ments of the density matrix, yielding the diagonal density 
matrix 



dec 
rran 



(n\a)\H,, 



2, _ p -M 2 H^ 



(5) 



The particle number in this mixed state follows the Pois- 
son distribution, and the average number of particles is 
n = \a\ 2 . The entropy content of the mixed state is given 
by 
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e " 2 — -(nhxn — n — Inn!), (6) 
n=o n ' 

where the superscript "cs" indicates that the result holds 
for a coherent state. With the help of Stirling's formula 
and the integral representation of the logarithm, 



In n = / — (e 

o s 



(7) 



the sum in © can be performed yielding an analytical 
result that is valid asymptotically for h S> 1 (actually, 
the approximation is excellent already for n « 1): 



~,(cs) 
^dec 



i(ln(2^) + l-i, + 



(8) 



It is not surprising that the entropy is proportional to 
In \ffi, because we have deleted all information about the 
relative signs of the amplitudes (a\n) by eliminating the 
off-diagonal elements of the density matrix. The num- 
ber of significantly contributing elements is given by the 
width, An = ^/fi, of the Poisson distribution. 

Let us mention that the energy for a single quantum 
oscillator in equilibrium at temperature T is given by 



ln(n + 1) + n In 1 + 



1 



(9) 



where n = [e^l T — 1) 1 is the average occupation num- 
ber. Asymptotically, for large n, one obtains S COl ~ 
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FIG. 1: Decoherence entropy Sdoc for a coherent state of a 
single field mode and equilibrium entropy S eq for the same 
average total energy as a function of the average occupation 
number n. 



^^dcc ' i- e - the thermal entropy becomes twice as large as 
the decoherence entropy. However, for small to moderate 
occupation numbers the ratio / S cq is close to unity 
and remains above 0.75 up to n = 10. Figure ^ shows 
the decoherence and equilibrium entropies as a function 
of the average occupation number n. It is evident that, 
for not too large values of n, the decoherence process 
generates a large fraction of the entropy that can be cre- 
ated, and any subsequent equilibration process adds only 
a small amount of entropy to it. Since decoherence is usu- 
ally a much faster process than thermal equilibration, our 
result implies that the fast entropy production observed 
in heavy ion collisions may be primarily due to decoher- 
ence of the initial state color fields. 

What does this imply for the quantum field theory, 
where the field is a system of (infinitely) many coupled 
oscillators? Assume that, after decoherence, the system 
can be described as a collection of N particles, given 
by some distribution function over single-particle states, 
which were generated by the decoherence of N cs coherent 
quantum states. Examples of such states are the internal 
wavefunctions of nucleons forming a large nucleus, or a 
quark with its comoving gluon cloud. Each coherent state 
contributes on average n — N/N cs partons. Then, after 
full equilibration, the thermal entropy is of the order of 
5th ~ N cs h = N, while for the decoherence entropy we 
get S'dcc ~ A^cgi ln(27m). The ratio of the two entropies 
is 



5th 



ln(27m) 
2n : 



(10) 



i. e. for large amplitude quantum states, which turn into 
many particles per coherent mode, the decoherence con- 
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tribution to the thermal entropy is small. On the other 
hand, if the individual occupation numbers are of order 
one, the contribution is sizable. This case applies to our 
problem of interest, the collision of two nuclei at high 
energy, as we will discuss below. 

For the coherent color fields in colliding nuclei, the 
average number of decohering gluons per transverse area 
has been given by 



dN _ C ¥ \n2Q 2 s Ay 
d 2 b TT 2 a> s 



(11) 



where Q s is the so-called saturation scale, Cf = 4/3 is 
the quadratic Casimir operator in the fundamental rep- 
resentation of SU(3), and Ay is the rapidity interval over 
which the color fields retain their coherence. The char- 
acteristic transverse area, over which the color fields are 
coherent, is ir/Q 2 , and one can argue that the longitudi- 
nal coherence length is of the order of Ay « l/a s ^3 • 
We thus obtain an average number of decohering partons 
per coherence region: 



C F In 2 



(12) 



For this value, our arguments presented above indicate 
that the entropy produced in the decoherence process is 
about half of the equilibrium entropy. The total entropy 
per unit rapidity produced by decoherence in a Au+Au 
collision at the Relativistic Heavy Ion Colider is 



dSdc< 
dy 



QjR 2 
2Ay 



(ln(27m) + 1) 



In 



2C F In 2 



1 



1500, (13) 



where we used the values [g Q 2 S ~ 2 GeV 2 , R = 7 fm, 
and a s s» 0.3. This value accounts for about half of the 
entropy measured in the final hadron distribution. 

Entropy production by decoherence of classical color 
fields was also discussed recently by A. Mueller on the 
basis of somewhat different arguments 01 > wn0 obtained 
a similar expression for the total generated entropy: 



2 D 2 



^dec ~ CsQ" s R 



(14) 



where cs is a nontrivial factor numerical of order one. 
It is the factor Q 2 S R 2 , which makes entropy generation 
by decoherence a large effect in his treatment, as well. 
Our results for the numerical factor differ in their depen- 
dence on a s , but not in their order of magnitude. An 
unambiguous determination of the factor cs will require 
a determination of the (de-) coherence length for the clas- 
sical fields, which exist in the nuclei before the collision. 

We finally discuss the case of particle production by a 
fast moving electric charge, such as a large nucleus. The 
occupation number of the various field modes forms the 
basis of the Weizsacker- Williams (WW) approximation 



to interactions of charged particles and photons at high 
energies 0]. For a Coulomb charge Ze moving close to 
the speed of light with Lorentz factor 7, the spectrum of 
equivalent photons is given by 



dn 2Z 2 a , / 7 
— = In I — — 

dw txuj \ ujR 



(15) 



where R is the intrinsic size of the charge and u < w max = 
7/i?. In order to determine the occupation number for a 
photon energy u), we need to know the coherence inter- 
val Aw, corresponding to the inverse longitudinal length 
scale of the process which causes the destruction of the 
coherent field. Characteristically, Aw grows with u>. For 
our analysis we assume that Aw ~ clo with a constant 
parameter e << 1 and discretize the integral over uj 



by setting w 7 - 



^e ej with integer j = 1, . . . , J and 



J = e _1 ln(w m i n /ti; max ). The entropy generated by deco- 
herence of the initial field configuration is then obtained 
as 



I 1 



i=i 
J 



4Z 2 aeln [ e~ tJ 



J2lH±Z 2 ae 2 (J 



i=i 

iln^.lnl 42T-,., Mu 

2e W m in 



(16) 



Thus the produced entropy depends crucially on the 
value of e, and thus on the process leading to decoherence. 
In the case of coherent color fields, a more quantitative 
understanding of the coherence length of the fields before 
and after the collision is desirable. 

In conclusion, we have shown that the decoherence of a 
quasiclassical state generates a significant amount of en- 
tropy. If the average occupation number of each coherent 
domain of the initial state is not much larger than one, 
the entropy released by the decoherence process is a siz- 
able fraction (e. g. one-half) of the entropy attained after 
thermodynamic equilibration. If this reasoning is applied 
to relativistic heavy ion collisions in the framework of the 
color glass condensate model, a significant fraction of the 
measured entropy of the final state can be produced on 
the time scale of decoherence, Td cc ~ 1/Qs < 0.2 fm/c. 
This may explain the rapid transition to a state that be- 
haves approximately like an equilibrated QCD plasma. 
The precise value of the entropy generated by decoher- 
ence Ijl3|l depends sensitively on the coherence length of 
the entropy creating process. 
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